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Abstract — Performance evaluation of the 802.11 MAC protocol 
is classically based on the decoupling assumption, which hypoth- 
esizes that the backoff processes at different nodes are indepen- 
dent. A necessary condition for the validity of this approach in 
the asymptotic sense (when the number of wireless nodes tends 
to infinity) is the existence and uniqueness of a solution to a fixed 
point equation. However, it was also recently pointed out that this 
condition is not sufficient; in contrast, a necessary and sufficient 
condition is a global stability property of the associated ordinary 
differential equation. Such a property was established only for a 
specific case, namely for a homogeneous system (all nodes have 
the same parameters) and when the number of backoff stages is 
either two or infinite and with other restrictive conditions. In this 
paper, we give a simple condition that establishes the asymptotic 
validity of the decoupling assumption for the homogeneous case. 
We also discuss the heterogeneous and the differentiated service 
cases and formulate a new ordinary differential equation. It is 
shown that the uniqueness of a solution to the associated fixed 
point equation is not sufficient; we exhibit one case where the 
fixed point equation has a unique solution but the decoupling 
assumption is not valid in the asymptotic sense. 

Index Terms — Mean field theory, ordinary differential equa- 
tion, fixed point equation, 802.11, decoupling assumption. 



I. Introduction 

IN order to live up to the growing and insatiable hunger 
for higher wireless throughput of users, IEEE 802.1 In 
was ratified and released in October 2009, to which a lot 
of enterprises are reported to have migrated. The increased 
maximum bit rate of 802.1 In, 600Mbps, along with its easy 
deployability, suggests the potential use of an 802.1 In access 
point as an wireless router transacting a huge amount of data 
of many nodes. In this work, we focus on the performance 
evaluation of 802.11 under the many-node regime as the 
population size (the number of wireless nodes) TV tends to 
infinity. 

Most existing work on performance evaluation of the 802. 1 1 
MAC protocol 0, flU, £□], ED relies on the "decoupling 
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assumption" which was first adopted in the seminal work by 
Bianchi Q. Though having been defined in various ways, it 
essentially assumes that all the nodes in the same network 
experience the same time -invariant collision probability, which 
in turn amounts to the assumption that the backoff processes 
are independent This assumption is unavoidable primarily 
because the stationary distribution of the original Markov 
chain cannot be explicitly written due to the irreversibility of 
the chain 1111 even for small number of backoff stages, i.e., 3 
and 4, unless the population is very small. A similar point was 
stressed by P. R. Kumar in an interview with Science Watch 
Newsletter lfl2l : 

"A good analogy is in thermodynamics. Instead of trying to 
study the behavior of just three or four molecules and how they 
move around, you study the behavior of billions and trillions 
of molecules. . . . Similarly, we want to see what you can say 
about wireless networks in the aggregate." 

which suggests an analogy of the intractable small-scale 
problems in different areas. If we liken each wireless node to 
a particle in a physical system, which condition would suffice 
for every particle being absolutely decoupled from the rest? 

Once we assume that the decoupling assumption holds, the 
analysis of the 802.11 MAC protocol leads to a. fixed point 
equation (FPE) ifTTI . also called Bianchi's formula. Kumar 
et al. IfTTI revisited the FPE and made several remarkable 
observations, advancing the state of the art to more sys- 
tematic models and paving the way for more comprehensive 
understanding of 802.11. Above all, one of the key findings 
of H2, already adopted in the field HH, Q3, is that the 
full interference model, also called the single-cell model fTTI 
and is the main focus of our work, leads to the backoff 
synchrony property lfl6l which implies the backoff process 
can be completely separated and analyzed solely through the 
FPE technique. 

It is unmistakable that the decoupling assumption is not 
correct for finite population (N < oo) due to the dependency 
between nodes, yet there is no theoretical support to whether 
this assumption is correct or not in the asymptotic sense as the 
population N goes to infinity. It is recently pointed out by the 
work by Benai'm and Le Boudec (H Section 8.2], using a mean 
field convergence method, that the uniqueness of a solution to 
the FPE does not necessarily lead to the asymptotic validity 
of the decoupling assumption in the stationary regime. That 
being said, the decoupling assumption along with its FPE has 

The meaning of "to decouple" in the literature as well as in our work is 
an abuse of terminology, in the sense that it has implied not only 'to decouple 
nodes' (independence) but also 'to have a time-invariant collision probability'. 
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not been theoretically validated yet in the asymptotic sense, 
though regarded as a plausible and intuitive approximation 
especially when the number of nodes is large (N ^> 1). The 
main purpose of this paper is to provide an answer to the 
following challenging question with appropriate mathematical 
formalism. 

"Under which conditions is the FPE asymptotically valid?" 

To put it another way, we revisit Bianchi's assumption upon 
mean field theoretic results (TJ, (6), |fl9l which state that, as 
N tends to infinity, a scaled version of the original Markov 
chain model of the backoff process in 802.11 MAC protocol 
converges to a nonlinear ordinary differential equation (ODE) 
so that the asymptotic validity of the FPE is rendered down 
into the stability of this ODE. Denoting by pf. the attempt 
probability of each wireless node at each time-slot in backoff 
stage k € {0, 1, • • • , K}, we assume in what follows that our 
mean field models are derived when K is finite and fixed while 
the number of nodes N goes to infinity. 

In connection with the mean field models, it is worth while 
to clarify why we use a specific intensity scaling regime, under 
which the activity of each node in backoff stage k is scaled 
as follows: 

p k := e(N) ■ q k (q k is a constant.) (1) 

where q k is called the scaled attempt rate throughout this 
paper. It is natural to assume that e(N) is vanishing, i.e., 
limjv^oo = 0. Otherwise, the collision probability 

between wireless nodes converges to one as N goes to 
infinity. More importantly, we have to use an appropriate 
form of intensity scaling e(N) in order to avoid exceptional 
cases. For example, if e(N) decreases faster than 1/N (e.g., 
e(N) = 1/N 2 ), it can be easily seen that the collision 
probability vanishes as N tends to infinity, irrespective of 
whichever backoff stage each node belongs to (we refer to 
Section ITl-El for a formal argument). In other words, each node 
is completely decoupled from the rest. On the other hand, if 
e(N) decreases slower than 1/N (e.g., e(N) = 1/yfN), the 
collision probability becomes one as N goes to infinity. That 
is to say, e(N) = 1/N is the only intensity scaling regime (up 
to a constant factor) that deserves to be analyzed. 

Under the intensity scaling regime e(N) = 1/N, Bordenave 
et al. in Theorem 5.4] studied the homogeneous case 
(all nodes have the same per-stage backoff probabilities) for 
the case when the number of backoff stages is infinite. They 
found the following sufficient condition for global stability of 
the ODE, hence for the asymptotic validity of the decoupling 
assumption: 

9o <ln2 and q k+1 = q k /2, Vfc > (BMP) 

where q k is the scaled attempt rate in (HJ for a node in 
backoff stage k. In this paper, we focus on the case where 
the total number of backoff stages K + 1 is finite, as this is 
true in practice and in Bianchi's formula. Sharma et al. |[T9l 
obtained a result for K = 1 and mentioned the difficulty to go 
beyond. A comprehensive summary of the literature and the 
outstanding questions raised therein has been recently made 
by Duffy 0. 



In the course of searching for the main result, we discover an 
intriguing fact that not only (i) the mono tonicity ( dMONOt in 
Sec. HIl) but also (ii) the mild intensity of scaled attempt rates 
(( IMINTI i in Sec. implies the uniqueness of a solution to 
the FPE, which is naturally a necessary condition for stability. 
Moreover, the latter ( IMINTb is proven to guarantee the global 
stability of the ODE. This finding greatly simplifies the story 
because simply that the attempt rate is upper-bounded by 
the reciprocal of the population suffices for each node being 
completely decoupled from the rest, namely q k < 1 for all k. 
Moreover, it is shown that, for the familiar parameter setting 
Qk = qo/m k where m > 1, the condition ( IMINTI ) suffices for 
maximizing the aggregate throughput of the network, hence it 
is a practical condition. 

In order to offer various services to higher priority users 
with additional performance requirements, 802. lie standard 
introduced the enhanced distributed channel access (EDCA) 
functionality that has two mechanisms to differentiate the per- 
class settings of (i) the contention window (CW) and (ii) the 
idle time after each transmission. Since the former, CW dif- 
ferentiation, necessarily implies there are two or more classes, 
we call the corresponding system heterogeneous. The latter, 
called AIFS differentiation, imposes an additional complexity 
on the Markov chain analysis because whether the users of 
a class may attempt transmission at each time-slot depends 
on the type of the current time-slot, which again depends on 
the activity of the users in the previous time-slot. This mutual 
interaction of the two evolutions substantially complicates the 
analysis. As of now, there is no ODE in the literature which 
models AIFS differentiation using an appropriate formalism. 

To tackle this problem, it is of importance to observe that the 
stage evolution of all nodes (or stage density) is much slower 
than the evolution of the type of time-slots under the AIFS 
differentiation. Thus the former can be taken to be constant 
by the latter. An application of mean field theoretic result HI 
Theorems 1 & 2], formalized based on the same observation, 
yields an extended ODE model of the backoff processes 
in EDCA-enabled 802.11 networks. We also formulate an 
extended FPE on the basis of this ODE, which is satisfied 
by the equilibrium points of the ODE. It is remarkable that 
this FPE coincides with that proposed in ifTTl Section VI]. 

The versatility of the ODE model is demonstrated by inves- 
tigating some selected counterexamples. In the first example, 
we consider a homogeneous system where all nodes use the 
same parameters and show that the system is bistable in 
that the backoff process, after whirling closely around an 
equilibrium for a very long time, suddenly jumps into another 
equilibrium, and vice versa. The FPE model is only capable of 
identifying three equilibrium points as its solutions, whereas 
the ODE model is further capable of classifying the two of 
them into locally stable points and the other into unstable 
point, accurately reflecting the multistability. The trajectories 
of the ODE constitute a separatrix which divides the initial 
condition space into two regions. We also consider a het- 
erogenous system where the set of nodes are divided into two 
classes. A delicate determination of the parameters renders the 
system oscillatory such that all trajectories converge to a stable 
limit cycle formed around an unstable unique equilibrium 
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point where the limit cycle is as determined by the extended 
ODE. This example also serves as an illustration of the fact 
that there may be a unique solution to the fixed point equation 
whereas the decoupling assumption does not hold in the 
asymptotic sense. 

This paper discovers a fundamental condition for perfect 
decoupling of wireless nodes under the asymptotic regime 
which is essential for the design of 802.11 in a predictable 
manner and an extended fixed point equation of greater 
generality which is founded upon mean field theory, which 
constitutes the contributions of the paper. We also stress that 
the stability condition established in this work for the first 
time has been tantalizing other researchers as well, e.g., |fl9l 
Appendix B]. 

The rest of the paper is organized as follows. In Section 
IIT1 we present a brief overview of recent advances in mean 
field theory and introduce the associated ordinary differential 
equation thereof. In Section [III] , we prove a global stability 
condition of the ODE, which is in turn shown to be capable 
of optimizing the throughput. In Section IIVI we elaborate 
on another complexities arising from EDCA and derive its 
corresponding ODE model. Some counterexamples in Section 
[V] illustrate the utility of the ODE models. Concluding remarks 
and an outstanding problem are given in Section [VI] 

II. Mean Field Technique Revisited 

The backoff process in 802.11 is governed by a few rules 
if the duration of per-stage backoff is taken to be exponential: 
(i) every node in backoff stage k attempts transmission with 
probability pk for every time-slot; (ii) if it succeeds, k changes 
to 0; (iii) otherwise, k changes to (k + 1) mod (K + 1) 
where K is the index of the highest backoff stage. Markov 
chain models, which have been widely used in describing 
complex systems including 802.11, however, very often lead 
to excessive complications as discussed in Section U In this 
section, we present a surrogate tool for the analysis, mean 
field theory. It is noteworthy that the rules used in 802.11, 
i.e., (i)-(iii), closely resemble the mean field equations laid 
out below. 

Note that our analytical model of 802.11 MAC protocol, 
from which we derive a set of ordinary differential equations 
based on the mean field argument, is different from the original 
one. Thus we first briefly describe the original operation of 
802.11 MAC in Section III-AI and explore the differences 
between our model and the real 802.11 MAC protocol in 
Section Illll 

A. Basic Operation of DCF Mode 

Time is slotted. Since our analysis is mainly focused on the 
backoff procedure of 802.11 distributed coordination function 
(DCF), we call the standardized time interval in the backoff 
procedure of the 802.11 standard time-slo^ for brevity. The 
durations of frames, packets, and inter-frame spaces used in 
the other procedures are generally different from that of a 
time-slot. 

2 This is equivalent to slot in the work by Kumar et at 1111 {e.g., 20fis in 
IEEE 802.11b). 



Each node follows the randomized access procedure of 
802.11 DCF. To begin with, each node generates a backoff 
value if it has a data packet to send. Since the backoff 
procedure of each node is controlled by inter-frame spaces that 
fill in spaces between frames and packets, we introduce them 
here to help to understand the basic operation of DCF mode. 
Two types of inter-frame spaces are used in 802.11 DCF, 
namely, Short Inter-Frame Space (SIFS) and Distributed Inter- 
Frame Space (DIFS). Each node freezes (stops) the countdown 
procedure as soon as the medium becomes busy. On the other 
hand, only when the medium is idle for the duration of a 
Distributed Inter-Frame Space (DIFS), a node may unfreeze 
(start) its countdown procedure of the backoff and decrements 
the backoff by one per every time-slot. If the backoff reaches 
zero, the sender transmits an RTS (ready to send) frame, 
followed by a CTS (clear to send) from the receiver, a data 
packet from the sender and an ACK packet from the receiver 
if the CSMA/CA (carrier sense multiple access with collision 
avoidance) is implemented. Note that SIFS must be smaller 
than DIFS so that no node is allowed to interrupt a sequence 
of frames and packets which are spaced out SIFS apart. 

There exist K+l backoff stages whose indices belong to the 
set {0, 1, • • • , K} where K > 0. If a node has not attempted 
transmission for a data packet yet, the node is supposed to be 
in the initial backoff stage where the backoff value is drawn 
uniformly from {0, 1, • ■ ■ , 2bo— 1} (or {1, 2, ■ ■ ■ , 26o}). Here 
2bo is the contention window that serves as the initial value of 
a backoff countdown. If two or more wireless nodes finish their 
countdowns at the same time-slot, there occurs a collision be- 
tween RTS frames if the CSMA/CA is implemented, otherwise 
two or more data packets collide with each other. If there is a 
collision, each node who participated in the collision multiplies 
its contention window by the multiplicative factor m = 2. In 
other words, each node changes its backoff stage index k to 
k+l and adopts a new contention window 2bk+i = 2m k+1 bo. 
If k + 1 is greater than the index of the highest backoff stage 
number, K, the node steps back into the initial backoff stage 
and the contention window is set to 26rj. 

Let L and L c denote the total durations of a successful 
packet transmission and a collision, expressed in terms of 
backoff time-slot. Also the fixed overhead for each successful 
transmission is denoted by L a . Note that L, L c and L a do not 
need to be integer numbers but can be arbitrary positive real 
numbers. In 802.11 DCF, if CSMA/CA is used, L represents 
the time to transmit an RTS frame, a CTS frame, a data packet, 
and an ACK packet plus inter-frame spaces, i.e., SIFS and 
DIFS, where L Q is L minus the time to transmit a data packet. 
The duration L c , much smaller than L, is the time to transmit 
an RTS frame plus one DIFS. 

B. Differences between Our Model and 802.11 DCF Mode 

Our analysis is made tractable by a number of differences 
between our model used in this paper and the original opera- 
tion of 802.11 DCF mode. First of all, we take the duration of 
per-stage backoff to be exponential as we did at the beginning 
of Section [D] Secondly, the parameter set is fixed for each 
version of the standard whereas each parameter in our model 
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may be an arbitrary number. For example, in the IEEE 802. 1 lb 
standard, m = 2, K = 6(7 attempts per packet), and 2&o = 32 
are used. 

We also make a few assumptions to simplify the exposition. 

• Single-cell assumption: Most importantly, this work fo- 
cuses on the performance of single-cell 802.11 networks 
in which all 802.1 1 -compliant nodes are within such a 
distance from each other that a node can hear whatever 
the other nodes transmit. Since all nodes freeze their 
backoff countdown during channel activity, the total time 
spent in backoff countdowns up to any time is the same 
for all nodes. Therefore, it is sufficient to analyze the 
backoff process in order to investigate the performance 
of single-cell networks. This technique has been adopted 
in many works including El, |6), ifTTI . ifTTll . 

• Greedy-node assumption: Secondly, we only consider 
the case of greedy wireless nodes that persistently con- 
tend for the wireless medium. 

• Error-free channel assumption: Lastly, we assume that 
the wireless channels are error-free so that failed trans- 
missions are caused only by collisions between RTS 
frames (for the case of CSMA/CA) or data packets. 

C. Bianchi's Formula 

In performance analysis of 802.11, Bianchi's formula and 
its many variants are probably the most known |2), 0, iflOl . 
ED, CI, 113, ED, HI. Assuming that there are N nodes, 
Bianchi's formula can be written compactly in a more general 
fixed point equation (FPE) form: 



Z^fc=o Pfc 

7 =l-(l-p)"-i (3) 

where p and 7 respectively designate the average attempt 
probability and collision probability of every node at each 
time-slot. The attempt probability in backoff stage k is denoted 
by pk and defined as the inverse of the mean contention 
window, i.e., pk = l/(bk — 1/2). Note that, as long as the 
backoff stage k = follows backoff stage k = K for any 
attempts, the statistics like p and 7 are not affected by whether 
attempts in the highest backoff stage K are successful or not. 

The FPE model has been used as a de facto principal tool 
for the analysis of the 802.11 MAC Protocol. The weak point 
of the FPE model is that it cannot be concluded entirely from 
the form of FPE whether its solution (even if it is unique) 
might be a good first-order approximation of p and 7. 

Exactly under which condition the FPE holds is recently 
being rediscovered with rigorous mathematical arguments El, 
l6l , ED, called mean field independence. This fundamental 
approach was originally developed in the work by Sharma 
et al. Ifl9l where the first mean field analysis of the 802.11 
MAC protocol was performed. Later, Bordenave et al. 
provided a broader mean field framework which extends to 
multiple-cell networks (cf. single-cell assumption in Section 
II-B) and supports the notion of 'resource' . The particle inter- 
action model proposed in JT] overcomes some limitations and 



broadens applicability of the model proposed in f6). The three 
works El, 0, lfl9l have found that, as the number of particles 
goes to infinity, i.e., N —> 00, the stage distribution of every 
node evolves according to a set of K+l dimensional nonlinear 
ordinary differential equations (ODE) under an appropriate 
scaling of time. 

D. The Mean Field ODE model 

Let us dive into the details of our 802.11 model and how 
the Markov chain of the model converges to the associated 
ordinary differential equation. 

Model description: In our version of 802. 1 1 DCF mode under 
the assumption made in Section III-BI there are N wireless 
nodes evolving in a finite state space {0, ■ ■ ■ , K } at discrete 
time-slots t G {0, 1, • ■ ■ }. Denoting by X n (t) G {0, • ■ ■ , K} 
the backoff stage (the state) of node n € { 1 , • • ■ , N} at 
time-slot t, we collect the observations X n (t), for all n € 
{1, • • ■ ,N} and compute the relative frequencies, which is 
called the occupancy measure (or empirical measure). For- 
mally, the occupancy measure in backoff stage k at (discrete) 
time-slot t is defined as 

1 N 

n=l 

where lr.\ is the indicator function. Let A T be the transpose 
of a matrix A. It can be readily observed that the occupancy 
measure vector <&(t) := ($o(i) ••• &K(t)) T possesses the 
Markovian property because all nodes in the same backoff 
stage are exchangeable under the greedy-node assumption 
in Section III-BI Thus the system can be described by K- 
dimensional vector €>(£) rather than A^-dimensional vector 
X(t) := (Xi(t) ■■■ X N (t)) T though the nodes are not 
distinguishable any more. 

This Markov chain (discrete-time Markov Process) is in fact 
analogous to the special continuous-time Markov process, i.e., 
density dependent population process, that was used in the 
seminal work by Kurtz J5J Chapter 11]. Basically, the two 
works, El and 0, are nontrivial extensions of the result in 
15) to Markov chain version. 

Key scalings: Unlike the density dependent population pro- 
cess in 1131 . our Markov chain in (|4]i cannot converge to an 
ODE as N — > 00 because a Markov chain evolves at discrete 
time-slots t € {0, 1, • • • }. The ODE is derived by means of 
the following two key scalings. 

• Intensity scaling is to slow down the evolution of each 
node by a factor of e(N), such that each node in 
backoff stage k attempts transmission with probability 
p k = e(N) ■ q k . 

« Time acceleration is to accelerate the evolution of time- 
slots by l/e(N), such that a variable at t before this 
operation is translated into another variable at t ■ e(N). 

Both El an d Jll rely upon the concept of intensity scaling 
for mathematical tractability. The main purpose of using the 
intensity scaling e(N) is to make sure that the intensity, 
defined as the number of state (backoff stage) transitions per 
node per time-slot, vanishes, i.e., converges to as N — > 00. 
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In our context, the intensity is pk = e(N) ■ qk, and thus we 
require 

lim e(N) = 0. 

N->oo 

In Section III-EI the implications of intensity scaling will be 
explored in detail in conjunction with the collision probability 
and its physical meaning. 

Since the intensity in the above vanishes, the number of state 
transitions of all nodes per time-slot is order of N-e(N) which 
is dominated by N, That is, the expected change of <l?fc(t) over 
two consecutive time-slots is order of e(N) which tends to zero 
as N — > oo. However, if we accelerate the evolution of time- 
slots by l/e(N), the change of $fc(t) becomes order of one, 
and thus the time-slots get closer, hence the time continuity. 
The limit variables which we obtain by applying the time 
acceleration and the limit operation N — > oo are dubbed mean 
field limits (MFL) in this paper. 

To avoid notational confusion, we use capital Greek letters, 
<£>&(■) (or $(•)) and T, to denote the original variables and 
lower-case letters, <pk (or </>(■)) an d 7 to denote their MFLs. 
The ODE : The scaled version of the Markov chain con- 
verges to an ODE system as N — > oo. It is shown 
in 0| that, as N tends to infinity, *(i/e(iV)) = 
($ Q {t/e(N)) ■ ■ ■ $ K (t/e(N))) T converges in probability to 
4>{t) := (0o (t) ■■■ <pK {t)) T which is the solution of the ODE: 

ffi(t) = q(t) (1 - 7 (f)) - qoMt) + q K <f>K(th(t) (5) 

inflow from K 

which is the differential equation with respect to <fio(t) and 

^(t) = q k -x4>k-i{th{t) - q k Mt), (6) 
at 

which is the differential equation for k € {1, • • • ,K}. Note 
that we denote by 

K 

q{t):=J2QkMt) (7) 

fc=0 

the MFL of the average attempt rate and j(t) is the MFL 
of the collision probability to be defined very soon. It is 
important to note that the above system is degeneratJ^, because 
we also have a manifold relation 4>o(t) = 1 — X^feLi 4>k(t), 
which can be plugged into (|6]l to eliminate (©, whereupon we 
only need to consider the ii'-dimensional system @. We will 
use the reduced version © throughout this work to simplify 
the exposition. This system © will be called homogeneous 
because all nodes adopt the same parameter set qk and K. 

The differential equation © can be intuitively understood. 
For example, the first term and second term on the right-hand 
side in © are respectively the inflow caused by collisions in 
the (k — l)th backoff stage and the outflow caused by any 
attempts in the kth backoff stage. Note that the underbraced 
term in © was not considered in 0, and exists only in 
networks with finite backoff stages. 

Collision probability: The full derivation of the ODE is 
omitted due to the space limit and a detailed one can be found 

3 A degenerate system has a singular Jacobian matrix which means that its 
linearization cannot determine the local stability of the system. 



in the works by Sharma et al. |[T9l Section III] and Bordenave 
et al. Section 5]. However, we believe that the readers can 
grasp the main idea by looking into the derivation of the MFL 
of collision probability in the following. 

Pick a backoff stage k' € {0, ■ • • , K}. For any node in 
backoff stage k' , the collision probability of the node at time- 
slot t is given by 

K 

r^k^-l-il-eMq*)- 1 IJ (1 - ^ ' Vkf* kit) ■ (8) 

fe=0 

This is the probability that at least one other node attempts 
transmission at time-slot t. Here we can see that the term 
e(N)qk' vanishes as N — >• oo. Thus we can define the MFL 
of collision probability as follows: 

7 (t) := lim T(t/e(N),k'). 

N— ¥00 

We assume the following special intensity scaling regime 
throughout the rest part of this paper: 



It follows from the definition of exponential function 

lim (1 - x/N) N = exp(-x) (9) 

./V— >oo 

and the definition of q{t) in (|7]i that 

j(t) = 1 - e~ E *=o <?fc0fc( * ) = 1 - c~ m (10) 

Also, remark that T(t, k) depends on backoff stage k, whereas 
its MFL j(t) is common to all nodes. 

E. The Intensity Scaling Regime 

Here we expatiate upon our discussion on the intensity 
scaling e(N) in Section Q] Note that the expression (TTOb holds 
if and only if e(N) £ Q(l/N). When e(N) £ Q(l/N), 
e.g., e(N) = l/N 2 or e(N) = l/VN, we can see from the 
forms of (O and (0 that ^(t) becomes either zero or one 
because K is a finite constant. Summing up, if we consider 
e(-/V) ^ <d(l/N), the decoupling assumption is asymptotically 
valid, which is in line with what intuition tells us. 

In connection with the above discussion, the intensity scal- 
ing technique can be construed as an essential property that 
must be imposed upon all practical systems where particles (or 
nodes) share a common resource of fixed capacity 0. If e(N) 
decreases faster than l/N, e.g., e(N) = l/N 2 , the common 
resource is not used at all as population tends to infinity, hence 
no collision. On the other hand, if e(N) decreases slower than 
l/N, e.g., e(N) = 1/vN, the common resource is utterly 
squandered in attempting transmission as N tends to infinity, 
ending up with collisions all the time. Therefore, the physical 
meaning of e(N) = l/N is crystal clear. 

F. Equilibrium Points 

Equating the right-hand sides of © and to zero yields 
the following equilibrium points: 

(?k = —7 <Po, and O = — =^ 



6 



whereupon the backoff stage distribution of every node at the 
equilibrium can be computed as: 



7 



Ik 2^=0 qj 



1j 
*,K 



By plugging the manifold relation X)fc=o^ fc W = 1 mt0 me 
above, we can get the following fixed point equation in the 
stationary regime: 

~ _ Z^fe=o 7 



Z^fc=0 q k 



7 



(11) 



(12) 



Note that, under the intensity scaling regime, i.e., p k = qk/N 
and p = q/N, (O becomes 



lim 1 - (1 - p^' 1 = lim 1 - ( 1 - — 



N-l 



= 1- 0" 



which is identical to (flZb . That is, we do not need to distin- 
guish between (fT2l and (0 under the intensity scaling regime. 

III. Asymptotic Validation of Decoupling 

The theoretical limit of mean field analysis represented by 
© needs to be clearly understood. The nonlinear ODE model 
only implies that any node will be in backoff stage k with the 
common probability <j>k(t/N) under the asymptotic regime, 
A*" — > 00. The component ratio 4>(t), in general a time-varying 
solution of ©, is not guaranteed to be constant. Bordenave 
et al. Theorem 5.4] studied its global stability of the 
asymptotic case when K = 00, but the more practical case for 
finite K remains to be proved. The need of a proof for finite 
K is stressed in 0] pp.833] due to its practical implication. In 
line with this, by appealing to a Lyapunov function, Sharma 
et al. proved this for the case K = 1 where there are only 
two backoff stages fl9l Lemma 3], 

A. Main Results 

Before presenting the result for finite K in Theorem [T] we 
describe two different sufficient conditions for the uniqueness 
of the equilibrium. To simplify the exposition, we first define 
two conditions: 

qk is nonincreasing in k. (MONO) 

([TT)-([T2]i has a unique solution. (UNIQ) 

The following lemma holds as long as the right-hand side of 
(fT2l is increasing in q. That is, the lemma does not fully exploit 
the exponential form of (fT2l . It is remarkable that Lemma [T] 
was originally established by Kumar et al. flTTl Theorem 5.1]. 
We give a simpler alternative proof in Appendix [A] based on 
the method of mathematical induction. 

Lemma 1 (Monotonicity Implies Uniqueness) 
( IMONOb implies dUNIQJ . 

To present the second sufficient condition for the uniqueness 
of the equilibrium, we define another condition: 

q(t) < 1, Vi > 0. (13) 



As we are interested in global stability, we need to show that 
the solutions of (|6]i with any initial condition converge to the 
unique equilibrium. Recall that q(t) = J2k=o Qk^ki^t fr° m 
the form of which it is clear that ( fl~3T > holds for any initial 
condition <fi(0) if and only if 



q k < 1, Vfc. 



(MINT) 



We call this condition dMINTb which is an acronym for 
'Mild INTensity'. In a sense, we can interpret the intensity 
scaling regime p k = qk/N as a way of weakening the node 
activity. From this point of view, ( IMINTI ) implies we impose 
an additional constraint upon the node activity. 

Interestingly, the above upper bound on the scaled attempt 
rates also implies ( |UNIQ| i. This intermediate result is pre- 
sented here to shorten the proof of Theorem [U which is the 
final form of the result. 

Lemma 2 (Mild Intensity Implies Uniqueness) 
JMINTt implies ( jUNIQJ . 

Pro of: Putting q max := max fee { ... >K y q k , it is clear that 
( IMINTb is equivalent to g max < 1. First, we have 



spK k spK k 
- 2^fc=o7 < 2^fc=o7 _ , 

Multiplying the both sides of ( fTTT i by yields: 



qe 



-q _ l^k=0 7 



3 _ 1 l^k=0 I 

i k ' L 7 k ' y^°° ^k 

2^fe=0 q k 2^k=0 q k ^k=0 I 



(14) 



(15) 



where the last equality follows from (fT2l . i.e., e^^ = 1 — 7 = 
VSfco7 fe - The second factor of the last equation of (fT5l ) 
can be rearranged as 



K 



J2 ^(i ~ 7) = (1 + • • • + i K ) ~ (7 + ■ ■ ■ + 7 



K+l\ 



k=0 



1-7 



A'+l 



l-(l-e-«) 



q\K+l 



which is a decreasing function of q. As the first factor of 
the last equation of <TT~5T > is also a decreasing function of q, 
(15[ is decreasing in q. On the other hand, qc~^ is increasing 



in q G [0, 1] and the range of qe q is [0, c 1 ] 
decreases from qo at q = to 



Since (Q~5]) 



EtoCl-e" 1 ) 



-l\k 



l^k=0 



(l-c- 1 ) k 



at q = 1, it suffices to show that the above is less than or 
equal to c _1 . In the meantime, dMINTI i implies that the above 
is less than or equal to e _1 . Therefore, (fTTT i and ( TTZb have a 
unique solution. ■ 

Unlike Lemma Q] the forms of both (fTTT ) and ( TTZb are fully 
exploited for the proof of Lemma 12 Specifically, the fact that 
(7(1 — 7) is an increasing function in q over the interval [0, 1] 
is used for the proof of Lemma |2] 

So far we have shown that there are two sufficient con- 
ditions, dMONOl i and (IMINTK for the uniqueness of the 
equilibrium, ( |UNIQ[ ), which is naturally a necessary condition 
for the global stability. We now show that one of them implies 
the global stability in the following theorem, which also 
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(MOISO) 
\ 



Fig. 1. Logical relations between conditions. 



completes the logical relations between dMONOK ( |UNIQ| >, 
dMINTI ), and the global stability, as shown in the Venn diagram 
in Fig. Q] Since it is not yet clear if there exists any case when 
dMQNQI i holds at the same time as the associated ODE is 
unstable, a part of the set ( IMONOb is depicted by the dashed 
line in Fig. Q] 

It is remarkable that Lemma [2] is now rendered obsolete 
by the following theorem because the global stability of (0 
automatically implies ( |UNIQ[ ), as clearly depicted in Fig. Q] 

Theorem 1 (Stability Condition) 
( IMINTI ) implies the global stability of ©. 

Proof: Because q(t) is bounded, there exist q l and a 
sequence {r;} such that 



liminf q(t) 



lim q{n) = q 



Since <p(t) is a probability measure on a finite sample space 
{0, ■ • • , A'}, 4>(t) is tight J4). Appealing to this, we can pick 
a convergent subsequence {ti} such that \im t t ^ oo 4>k(ti) = 
<pk(oo) exists. 

Defining v(t) = inf s > t q(s), we necessarily have v(t) < 
q(t), \/t > and lim^oo v(t) = q l . Consider the degenerate 
version of © which has one additional equation with respect 
to £r(t). By replacing q{t) with v{t), we get the following 
modified ODE: 

^(t) - v{t) C ~^ - q oVQ (t) + QKVKiW - e-"«), 



dt 



(t) = qk-i<p k -i(t)(l - e-"W) - qk<Pk(t). 



Since v{t) becomes a constant for t = oo, this ODE reduces 
to a linear ODE as t — >• oo whose coefficient matrix takes the 
following form: 



/ 



-go 
qoi 1 
o 







V 



9i7 








91 

i 



'72 










-QK J 



where we used 7 := lini4_ > . c>0 (l — e~ u ^ ') for notational sim- 
plicity. Applying Gershgorin's circle theorem to the transpose 
of this coefficient matrix shows that all the eigenvalues are 



negative hence that <p(t) converges as t 
should satisfy 



00. Thus <p(oo) 



q (p (oo) = q l e ql + qK<PK(ca) f 1 - c q '^j , 
q k (p k (oc) = g fe _i</J fe _i(cx>) - e -5 J , 



(16) 
(17) 



for k £ {1, • • • , A"} because lim t _ s . oc , y(t) = q'. Plugging (TTvT > 
into (fTol l yields 



q k ip k {co) = q\l 




-1 v 



(18) 



Suppose the initial condition <pfc(0) = </>fc(0), Vfc 6 
{0, • • • , A}. We have the following equations from the mod- 
ified ODE: 



<Po(t) 

+ J* e 9«( s -*) {j/( S )c-^ s ) + q K <PK(s) (1 - C-»^) } ds, 



(19) 



e" 9fc Vfe(0) 



(20) 



¥>*(*) 

+ /„* ^('-Ogfc-i^-iCs) (1 - e-'W) ds, 

where k £ {l,'" >-K"}- Fi rst we have v(t) < 1 from the 
assumption ( IMINTt . Since 1 — e~ x and xe~ x terms in the 
above equations are increasing functions when x € [0, 1] and 
vft) < q(t), it can be checked by plugging ( f20b into (fT9l A' 
times that <po(t) < 0o(i) and hence (pk(t) < <fik(t), Vt > 
and Vfc 6 {0, ■ • • , A}. That is, <f>k(f) is lower-bounded by 
Pk(*). 

From ([ToT i and the definition of the subsequence {U}, we 
have the following relation: 



EfcLo qk¥k(oo) = q l = Efc=o qk4>k{oo). 

where we recall ^fc(oo) was defined as limt^oo 4>k{ti) = 
<pk{oo). This result taken together with tp k (t) < <f>k(t) 
proves ipk(oo) = <f>k(oo), Vfc <G {0, • • • , A}, and therefore, 
Y^k=o fk{oo) — 1. Then it necessarily follows that q l should 
satisfy (fTTl i and ( fT2l which have a unique solution by Lemma 
12 This implies q l = q. 

Note that we can also prove <f = q in a similar way 
by defining <f and {i;} such that limsup^^ q(t) = g", 
limt^oo q{ti) = q u and limt^oo (pk(U) = 0fc(oo). This will 
show limt^oo = q. That is, there is only one limit point 
for q(t). 

Finally, we can pick a new sequence {1^} such that 



K 



lim inf 4>k(t) 

t— )-oo 



4> l k , lim fc (ri) = <tj 



for all fc € {0, ■ ■ • , A}. Using the fact lim t _ ! . 0O q(t) = q, it can 
be easily proven that lim^co 4>kif) = 4>k, Vfc € {0, • • • , A}, 
in a similar way. This establishes that (0, q, 7) is globally 
stable. ■ 

Remark 1 This result gives an answer to the question raised 
in Section U and justifies the FPE approach used in J2], 
IE), (iQl, d], CD, CSl, H71, (HI under a special scaling 
regime. That is, the decoupling assumption is validated in the 
asymptotic sense, as long as the scaled attempt rates are mild, 
i.e., ( IMINTt . 
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The result of (6] Theorem 5.4] implies that, for the case 
K = oo, a set of strong conditions is required for the global 
stability of the ODE; a monotonicity condition along with a 
condition on attempt rate in backoff stage k = 0, as shown 
in the condition ( IBMPb in Section U These strong conditions, 
designated also by ( IBMPI ) in Fig. Q] were proven to prevent 
wireless node to escape to infinite backoff stage. We can see 
that they correspond to a proper subset of the intersection of 
( IMINTI ) and ( IMONOb . As compared with Theorem 5.4], 
Theorem Q] is a stronger yet more practical argument due to 
finite K. 

As shown in Fig. Q] while the monotonicity ( IMONOb 
implies only the uniqueness (UNIQi which is not a decisive 
factor, ( IMINTI i implies both ( jUNIQ| i and the global stability, 
assuring the asymptotic validity of the decoupling assumption. 
It is still open whether ( IMONOt implies the global stability 
or not. 

Informally, the proof of Theorem [T] follows from the fact 
that the solution <p(t) cannot have more than one limit point. 
The key observation underlying its proof is that there exists 
a stable differential equation which becomes asymptotically 
linear as t — > oo at the same time as its solution (p(t) lower- 
bounds 4>(t) such that <fi(t) is squeezed into <p as t — > oo. 

It is an intriguing fact that the above theorem may be 
restated in terms of 7(i) rather than q^, hence an alternative 
interpretation of the theorem: the ODE is globally stable if the 
collision probability -f(t) < 1 — c _1 = 0.632 for any initial 
condition 0(0) = (<j>o(0) ■■■ 4>k(0)) T - This interpretation 
means that if the collision probability is small enough, then the 
decoupling assumption is asymptotically valid, which appears 
to be in best agreement with our intuition. 

As was mentioned at the beginning of this paper, there is 
only one intensity scaling regime e(N) = 1 /N which deserves 
to be analyzed because, under this regime, it is not clear 
whether the collision probability would converge to a unique 
equilibrium point and would stay around there forever. Though 
we have also shown that ( IMINTI i is a sufficient condition 
for the asymptotic validity of the decoupling assumption, one 
may ask in return whether there exist any example where 
( IMINTI ) does not hold and the decoupling assumption is not 
asymptotically valid. Yes, there is. We will show in Section 
IV-BI that the collision probability may oscillate between two 
values as time goes if ( IMINTI i is violated. 

B. Achievable Throughput 

Recall that L and L c denote the durations of a successful 
packet transmission and a collision, expressed in terms of 
backoff time-slot. The fixed overhead for each successful 
transmission is denoted by L a . In what follows, we make a 
mild assumption that L c > 1 which means that the duration 
of a collision is no less than that of a single backoff time-slot. 
As we explained in Section IH-AI L c is an RTS frame plus a 
DIFS, both of which is larger than a backoff time-slot in all 
versions of IEEE 802.11 MAC, regardless of the usage of the 
RTS/CTS mechanism. 

Assuming that q(t) — > q as N tends to infinity, we can 
define the achievable throughput or alternatively the MFL of 



the aggregate throughput, as in J5J Section 5]: 



Q(q) 



Pl(g) • (L + L ) + P (q) + P c (q) ■ L c 



(21) 



where Pi(g) := ge"«, P (g) := c"«, and P c (g) := 1 - 
Pi(g) — Po(q) are the MFLs of the probabilities at each time- 
slot that only one node attempts transmission, none of the 
users attempts transmission, and at least two users attempt 
transmissions, respectively. Derivations of these MFLs are 
similar to that of (O and thus omitted. 

Since (f2TT> holds on the condition that ^ is globally stable 
such that limt-^oo q(t) = q, we can use (f2TT> so long as ( IMINTI ) 
holds. Then the result of Theorem Q] poses another question: 

"Is there qk satisfying ( IMINTI i and maximizing ( ETl ) as well?" 

Dividing the denominator of ( 1211 ) by its nominator, we can 
see that maximizing (l2T~b is equivalent to minimizing 

1 

-{l-L c ) + -L c . 

g q 

Differentiating this expression shows that the global maximum 
of (fSTJ is at the solution of the following equation: 



1 



-l = (g-l)e« 



(22) 



whose left-hand side is monotonically decreasing in L c over 
the domain (0, oo) and whose right-hand side is monotonically 
increasing in q over the same domain. Also both sides have 
the same range, i.e., (— l,oo). This implies, for each value of 
L c E (0, oo), there exists a unique solution to (1221 , which is 
from now on denoted by q = q*. If L c = 1, the solution is 
q* = 1. Putting these facts together, we can see that, if L c > 1, 
there exists a solution q* < 1 to ( l22l which maximizes (fJTJ. 

It is more important that qk satisfies (IMINTI i because 
(IMINTI i is a sufficient condition (and the only sufficient one 
we know) for the throughout equation (l2T~b to hold. To this aim, 
we show here that there are infinitely many constructions qk 
what satisfy ( IMINTI i and ( IMONOl i and maximize d2TJ at the 
same time. For q^ = qo/m k and q = q*, plugging (fT2l into 
(fTTT i yields: 



<1_ 



E 



K 

k=0 



1 



Ef= (i 



(23) 



The right-hand side of ( f23l is decreasing in m g (0, oo). For 
given optimal solution q*, one can use d23l to find go and m 
which satisfy (IMINTI i and maximize OTb at the same time. 
For instance, in order to obtain nonincreasing qk, one can 
simply set go = 1 an d compute m from d23l where m > 1 is 
warranted because the left-hand side of d23l is no greater than 
1 and the right-hand side of d23l decreases from 1 at m = 1 
to at m = oo. 

To sum up, for every go £ [g*, 1] where g* is the solution 
to d22l . the construction g^ = qo/(m*) h , where m* is the 
solution to d23l ). maximizes the aggregate throughput (|2~T1) as 
well as guarantees the global stability of ©, 



9 



IV. Mean Field with Service Differentiation 

So far the discussion has centered on the homogeneous 
system where all nodes have the same parameter set. Now 
we turn to the heterogeneous case arising from the service 
differentiation mechanisms defined in 802.1 le standard. In 
addition, a special kind of coupling caused by one of the 
mechanisms necessitates formulating a new ODE model. 

A. Prioritization Mechanisms 

Although three prioritization mechanisms are provided by 
enhanced distributed channel access (EDCA) functionality, 
one of which, called transmission opportunity (TXOP) Q, 
exerts its influence only on time-slots when all nodes are 
freezed (See Section Ill-Al l, hence no need for making an 
analysis of it. The other two mechanisms are to differentiate 
per-class settings of 

• contention window (CW), 

• arbitration interframe space (AIFS). 

The first mechanism, CW differentiation, in the present 
context amounts to per-class setting of qo and K, on the 
assumption that qf. = qo/2 k for k <G {0, • • • ,K}. We 
extend this feature by allowing per-class setting of K and 
qk for any k £ {0, • • • , K} for the sake of generality and 
notational aesthetics. Since CW differentiation implies that 
there are two or more classes, the corresponding system will 
be called heterogeneous, whether the following differentiation 
is enabled or not. 

The second, called AIFS differentiation, is to offer a soft 
non-preemptive prioritization to a certain class by holding 
back other classes from attempting transmissions for a few 
time-slots. This prioritization is effectuated by idling nodes 
for different durations, i.e., AIFS, after every transmission. In 
other words, AIFS differentiation reserves a few time-slots for 
high-priority classes. 

The analysis here is presented for the case where there 
are two classes, i.e., Class H (high) and Class L (low), only 
to simplify the exposition, but can be extended to arbitrary 
number of classes. Let us call the time-slots reserved for Class 
H reserved slots, which will correspond to the superscript R. 
We call the remaining slots following reserved slots common 
slots, corresponding to the superscript C. Note that both Class 
H and Class L users can access the channel during common 
slots, whereas the backoff procedures of Class L users are 
suspended during reserved slots. The per-class parameters and 
occupancy measures are denoted by q^, q]:, K H , K h , 
and $L(t). 

There are two kinds of couplings caused by the above- 
mentioned prioritization mechanisms. 

• Inter-class coupling: As compared with the analysis 
carried out in Section III-CI where the stage evolution 
of nodes depends only on their own stage density, i.e., 
the occupancy measure $(<) = ($o(t) ■■■ $/^(i)) T , 
the performance analysis of 802.11 in the presence 
of CW differentiation is complicated by the very fact 
that two-class users mutually interact with each other 
through * H (t) := •■• $^„(t)) T and * L (t) := 



($o(t) • • ■ $^. L (i)) T . Fortunately, it turns out not very 
difficult to incorporate this complication into the ODE 
model in the previous section because we are simply 
dealing with two evolutions of the same kind. 
• Coupling between two kinds of evolutions: However, 
when it comes to AIFS differentiation, the issue is in- 
volved by the fact that the stage distribution of nodes 
in the previous time-slot affects the type of the current 
time-slot, and besides, the type of the current time-slot 
also affects the stage distribution of nodes in the next 
time-slot. That is, there are now two different kinds 
of evolutions, stage evolution of nodes and slot type 
evolution of time-slots, the latter of which adds a new 
type of state variable to the Markov chain model in 
J2j. An interesting point to note is that Sharma et al. 
|fl9l Section IV] in a similar context also reckoned this 
difficulty though they have not solved it. 

B. Markov Model for the Evolution of Slot Type 

To avoid notational confusion, we use only the original 
occupancy measures in discrete time, i.e., 5>j?(t) and <S?t(t) in 
this subsection. The MFLs of these variables will be defined 
in the next subsection. We first divide the population into two 
classes such that 

N 



N + N = N, <j n :=—, <t l := — 



N L 
~N 



Without loss of generality, the sets of nodes of Class H and 
Class L are denoted by N H := {1, • ■ ■ , N u } and N L := {N H + 
1, ■ • ■ , N}. Thus we define the occupancy measures as 



ngN" n£N L 



{X n (i)=k} 



so that we have a 11 = ££!„**(*) and a L = ZCo*W)- 
Since there is no inter-class transition of users, cr H and cr L are 
constant and satisfy the relation <r H + er L = 1. In this setting, 
the probability that one or more nodes attempt transmission at 
time-slot t of slot type R or C is as follows: 



K 



r R (t):=i- [ 1 



k=0 



Ik 

N 



fe=0 



rc (i) : =1 -n(i-£ n i-% 



fe=0 



where we intentionally use the letter V which is the same as 
the collision probability in ^ because in the mean field limit 
the additional term q^ /N in ([H) vanishes, and thus the MFLs 
of the above equations and (© are much alike. 

From the viewpoint of an individual node, we can describe 
AIFS differentiation by only three rules: (i) after any trans- 
mission attempt which is either successful or a failure, AIFS 
procedure is initialized, i.e., a counter value is set to zero; (ii) if 
the current time-slot is idle, the counter value is incremented 
by one; (iii) if the counter value reaches its per-class AIFS 
value, the node may attempt transmission with its per-stage 
probabilities, i.e., q^- m/N and q 1 ^- r t )/N, 
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i-r E (t) i_r R (t) i-r E (t) i-r E w 



l-r^t) 



r"(t) 




" ^ Common Slots^ 
" ""PI . 1^ ""PI 



Fig. 2. Evolution of slot type follows a nonhomogeneous Markov chain. 



Denoting the difference of the two per-class AIFS values by 
A > 0, we can see that the transition structure based on the 
aforementioned rules are illustrated by the nonhomogeneous 
Markov chain in Fig. [2] where we used the non-idle proba- 
bilities, i.e., r R (i) and T R (i), and the idle probabilities, i.e., 
1 — T R (<) and 1 — r c (i), as well. Here in Fig. preserved 
time-slots and common time-slots are respectively denoted by 
the notations 'Slot l'-'Slot A' and 'Slot A+'. Note that A+ 
means that, after any A or more consecutive idle backoff time- 
slots, the corresponding slot-type must be C. It should be clear 
in Fig. |2] that not only slot-type but also the backoff stages of 
nodes, i.e., <l? H (i) and $ L (t), are also changing over time- 
slots, hence T R (t) and r c (t) are. 

The simplification of the analysis bases upon the following 
intuitive observation: 

"As population grows, the stage distribution (density) varies 
much slower than the type of time-slots." 

This observation follows essentially from the intensity scal- 
ing. Formally speaking, the occupancy measures, 3>j?(t) and 
5?fc(f), evolve at a rate of 0(1 /N) which ultimately vanishes 
as N — > oo, whereas the probability that the slot-type changes 
for each time-slot does not vanish and is strictly positive. 
Therefore, we can analyze the evolution of slot type as if 
the occupancy measures were constant. Solving the balance 
equations as if the Markov chain were homogeneous yields 
the following stationary distributions for each slot type: 



n c (i) 



Ef=o (i -!*(*))* 



{Ef=o 1 (i-rW} + ^ 



r R (t)) z 



rc(t) 



(i-r R (t)Y 



{Ef=o 1 (i-rW} + ^ 



rc(t) 



which satisfy n R (t) +U C (t) = 1. Note however that in general 
it is impossible to derive the stationary distribution of nonho- 
mogeneous Markov chains where the transition probabilities 
are time-varying. 

C. Extended ODE Model with Prioritization Mechanisms 

The MFLs of $£(t) and <$>\(t) are denoted by <f%(t) and 
4>\(t) as in Section IH-DI After manipulation akin to ( fTOl i. 
we can show that the MFLs of collision probability for the 
different types of time-slots, R and C, take the forms 
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R (t) = 1 



= i - e-? H W 



and 7 C (t) = 1 - c 



-q H W-q L (t) 



where the MFLs of the per-class average attempt rates are 
defined as 



fit) 



fe=0 



A 



and q h (t):=J2<lk<f>W)- 



k=0 



Let g R (t) and g c (t) denote the derivatives of (<p H (t) , 4> h (t)) 
under the conditions that all time-slots are of slot type R or C, 
respectively. By means of the informal observation in Section 
IIV-BI Benai'm and Le Boudec have formally established a few 
results (TJ, by appealing to which, we now derive the final 
extended ODE. To this aim, as in Section III-DI assume that 
we eliminate manifolds by using the equations <$$(t) = er H — 

E*Li^fc(*) and <#>(*) = aL - Ef=i and hence we 

consider (K H + K h ) -dimensional ODE. From the fact that at 
a time-slot of slot-type R, only Class H users are allowed to 
attempt transmission, whereas at a time-slot of slot-type C, all 
users are allowed to do so, we have 



9 R (t) = 



9 C (t) 



( 4h_ 1 ^h_ 1 W7 R W-<? A -h^hW \ 

«(*)7 R (*) 






Interestingly, we can apply |fl] Theorems 1 & 20 to show 
that the resultant derivatives of (<fi H (t), (fi h (t)) becomes: 

s(^))^ R<i) -' R<,)+! ' C( ' wC( " (24) 

where ir R (t) and 7r c (t) take the following forms 



{£to 1 (i-7 I W} 

(l-7 R (t)) A 



(1-7 R W) / 
7 C M 



7 C (*) 



{Ef= 1 (i-7 R (t)r} 



(WW 

7 U W 



This result has the implication that we can derive d24l i as 
a linear combination of the MFL vectors, i.e., g R (t) and 
g c (t), with the coefficients, i.e., n R (t) and n c (t), which can 
be defined as 

ir R (t) := lim IL R (Nt), ir c (t) := lim if(Nt) 

N— >oo JV— >oo 

where IT R (-) and II C (-) are the stationary distributions we com- 
puted from Fig. |2] in Section HV-Bl as if the nonhomogeneous 
Markov chain were homogeneous. 

4 The corresponding assumptions can be easily checked. 
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Finally, after some manipulation of (l24l i. we have the 
following enhanced ordinary differential equation: 



(25) 



dt 



■(t) = n c (t) {gti^-i (th c (t) - , (26) 



where d25l l and d26l i respectively hold for k £ {1,- • • ,-ftT H } 
and fc e {1) • " ) ^ L }- Here we use the following shorthand 
notation: 

^(t) = n K (t)^(t)+n c (th c (t) 

whose form is obvious from (l24l >. 

In the stationary regime, we can get the following fixed 
point equation: 



l^k = q« 



l^k=0 q L 



y H = ( ! _ e -5- 



7T C ll-a-r-T 



7 L = 1 - e" 



(27) 

(28) 

(29) 
(30) 



Remark 2 It is remarkable that the extended ODE model laid 
out in d25l l and ( T26l l encompasses the homogeneous system in 
Section HIl and the heterogeneous system in Section lTVl as well, 
which has the two prioritization functionalities. 

For instance, if A = 00, we have ir R (t) = 1 and the ODE 
model reduces to the homogeneous system ©. On the other 
hand, if A = 0, we have n c (t) = 1 and the ODE model 
reduces to a purely heterogeneous system, implying that the 
AIFS differentiation is disabled. 

What is the most surprising is that the FPE (|27|i-(|3"0l 
coincides with that proposed in ifTTl Section VI], which was 
derived rather intuitively. 

In the following, we introduce three new conditions akin 
to those in Section III-DI by adopting which we present two 
lemmas. 

and q k are nonincreasing in k, (31) 
(f27T>-(T30l> has a unique solution, (32) 
ql < 1 and q\ < 1, Vfc. (33) 

Lemma 3 (Monotonicity Implies Uniqueness) 

(ED implies d32l 

Lemma 4 (Mild Intensity Implies Uniqueness) 
< f33T > implies d32l 

Proofs of the above two lemmas are in Appendix. It is 
of importance to note that these lemmas are of even greater 
generality because they hold for all A > and A = 00 as 
well, implying that Lemmas Q] and [2] respectively correspond 
to the special cases of Lemmas [3] and [4] i.e., the case A = 00. 



We are not able to prove the equivalent of Theorem [T] for 
the case where there are more than one class due to the inter- 
class coupling arising from CW differentiation. This coupling 
makes it technically challenging to find a stable ODE, which 
would bound the solution of the ODE as in the proof of 
Theorem Q] In the meantime, the other coupling induced by 
AIFS differentiation does not seem to cause a major technical 
difficulty. As of now, we have to be content with having stated 
the problem precisely with its inherent technical difficulty. 

V. Selected Counterexamples 

Before proceeding to selected examples, the gap between 
the ODE model and the backoff processes in 802.11 must 
be bridged. This gap emerged right on applying the intensity 
scaling in Section Hl-DI The scaling relation pk = qk/N sug- 
gests to us that replacing qk by Npk should yield a reasonable 
approximation if pk is small. Plugging this approximation and 
removing the time acceleration from ©, we have 



dt 



(t) = pk-i(/)k-i(t)j(t) -pk<f>k(t) 



(34) 



where 7 (t) := 1 - e~ N m and p(t) ££=oPfc&(*)- 

We compare the trajectories of the ODE model d34b with 
the simulation results of the corresponding Discrete Time 
Markov Chains (DTMC). Since all Markov chains used for the 
simulation are ergodic: statistically speaking, the two system 
forget their initial states after evolving for a long time, we 
run each simulation for as long as possible. To this aim, we 
have run each of the DTMC simulations for 120,000,000 
backoff time-slots. It is remarkable that we must run the 
DTMC simulations for a long duration because the DTMC 
in Section IV-AI exhibits a very unique phenomenon called 
bistability {e.g., see Fig. |3(b)[ ). 

To obtain the short-term average statistics, the entire du- 
ration of each simulation is divided into disjoint intervals of 
2, 000 time-slots and each short-term average data point was 
calculated over one of the disjoint intervals. We assume that 
all wireless nodes are in backoff stage at the initial time-slot. 

A. Example 1: Multistability 

Consider the homogeneous system (l34l . Plugging (j2]i into 
© yields: 



/( 7 ) :=l-exp(-iV ^ ^ 



-7 = 



(35) 



which is a function of only 7. Consider the following mul- 
tistability example where there are N = 1200 nodes and 
K + 1 = 13 backoff stages. The attempt probability at each 
backoff stage pk is 



1 



1 



3200' 160' 160' ' 160 



{po,Pi, ■ ■ ■ ,pi 2 ) 



where m = 6/5 = 1.2. The roots of (1351 1 can be computed 
from Fig. |3(a)| as 

(71,72,73) = (0.540,0.828,0.952). 
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(a) f(~f) versus 7 




Backoff time-slots x 10 6 

(b) Short-term average collision probability vs. backoff time-slots 
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0.78 0.8 0.82 0.84 0.86 0.88 0.9 0.92 0.94 0.96 0.98 

Short-term occupancy measure of stage 

(c) Short-term average occupancy measure in backoff stages, 
(<7>o(*),<M*)). (Mt),<l>2(i)) and W>0 (*), <t>Z (*))• Stars (*) and 
circles (o): mean field limits; dots: DTMC simulation. 

Fig. 3. Bistabilty Example: There are three solutions to the fixed point 
equation, two of which (71 and 73) are stable and the other one (72) is 
unstable. Short-term average statistics measured for each 2000 backoff time- 
slots suggest bistability. 



The instantaneous collision probability for each 2000 back- 
off time-slots is shown in Fig. |3(b)[ which tends to concentrate 
around 71 = 0.540 and 73 = 0.952. Note that the average 
collision probability for the entire duration of the simulation 
is 0.832 that is neither 71 nor 73. Recall that <f>k(t) denotes 
the fraction of nodes in backoff stage k. Fig. |3(c)| shows the 



short-term average of the fraction of nodes in backoff stage 
k € {1,2,3} versus that in backoff stage 0. From the top 
to the bottom, the short-term occupancy measures of stages 
1 — 3 are shown in order, where the two kinds of markers, i.e., 
circle (o) and star (*), stand for the occupancy measures at 
two equilibriums, 73 and 71, which are computed from ( l34l . 
The bistability of this system is precisely predicted from either 
two modes of behavior of (l34l or the eigenvalues of Jacobian 
matrices at the three equilibrium points. 



B. Example 2: Stable Oscillation 

We have managed to discover a rare example by delving 
into the heterogeneous system, without AIFS differentiation, 
i.e., A = 0, which in turn leads to 7r R = and tt c = 1. 
Suppose there are two classes H and L such that population 
of each class is N H = N h = 640. The numbers of backoff 
stages are assumed to be equal, i.e., K H + 1 = K h + 1 = 21. 
The attempt probability at each backoff stage is: 



,P2 ) 



1 



1 



2400' 480' 40' 



,19 



40 



/ 1 



1 1 



V 3840 64' 64' 



1 

' 64 



where m = 4/5. It is easy to verify that the corresponding 
fixed point equation takes the following form: 



/(7) 



1 n exp 

Xe{H,L} 



-N 



x l^k=o 7 



7 = 



which has the following unique solution as shown in Fig. |4(a)| 

7 H = 7 R = 7 C = 71 = 0.912. 

Since there is only one solution, one might be much inclined 
to hazard the conjecture by Bianchi et al. J2], ifTTl that the 
collision probability is approximately 71. However, there is a 
stable limit cycle around this equilibrium. In other words, the 
oscillation is stable, i.e., not transient but lasting forever. The 
event-average collision probability obtained through simula- 
tions is 0.869 which is less than 7 H or 7°. 

We can see from Fig. |4(b)| that, unlike the previous example, 
the trajectory of instantaneous collision probability forms 
almost periodic oscillation and does not tend to concentrate 
around the unique equilibrium 71. Though the oscillation is 
not deterministic but stochastic, it clearly persists as time goes 
to infinity. The period of the oscillation empirically can be 
computed from Fig. |4(b)| as between 19000 and 20000 time- 
slots. The oscillation and its period are exactly predicted from 
the trajectories of the ODE model (sold lines) as shown in Fig. 
|4(c)| The unstability of 71 can be decided by the eigenvalues 
of the corresponding Jacobian matrix. 

The decoupling assumption does not hold in the asymptotic 
sense; in contrast, nodes are coupled by the oscillations of 
the occupancy measure, an emerging property of the system 
dynamics. 
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(c) Short-term average occupancy measure in backoff stages, 
(<j>o(t), 4>i (t)) and (if>o(t), (f>n(t)). Solid lines and stars (*): mean 
field limits; dots: DTMC simulation. 

Fig. 4. Oscillation Example: There is a unique solution (71) to the fixed 
point equation but the decoupling assumption does not hold in the asymptotic 
sense. Short-term average statistics measured for each 2000 backoff time-slots 
suggest stable oscillation around the unique equilibrium. 



VI. Concluding Remarks with a Conjecture 

Since it is axiomatic that the fixed point equation (FPE), 
called Bianchi's formula, must have a unique solution in order 
to provide an approximation, there has been a speculation that 
the uniqueness of the solution might assure the validity of the 



FPE, which has been the main subject of previous approaches 
by Kumar et al. ifTTl . and Ramaiyan et al. ifTTl . One coun- 
terexample in our paper has shown that this speculation is not 
always true, putting another emphasis on asymptotic validation 
of the decoupling assumption which underlies the formula. 

Thanks to recent advances in mean field theory JT], (6) 
and also jTgj, we have analyzed the validity of the FPE by 
determining the stability of an ordinary differential equation 
(ODE). In the course of establishing stability, we obtained 
an illuminating insight that not only monotonicity but also 
mildness of scaled attempt rate guarantees the uniqueness of 
the equilibrium, which made the logical relations between 
them clear. Paradoxically, the mathematical formalism of mean 
field theory presented us a succinct stability condition (1MINTI ). 
whose main implication is as follows: to achieve perfect 
decoupling between nodes as population N grows, in addition 
to reducing the attempt probability at fcth backoff to qk/N, 
we need to further diminish the node activity such that the 
scaled attempt rate satisfies qk < 1. The existence of such an 
upper bound appears to be in best agreement with our usual 
intuition. We also have shown that there are infinitely many 
constructions of q^ which maximize the aggregate throughput 
as well as satisfy ( IMINTb and (iMONQl i. hence the condition 
(IMINTb is practical as well. 

Though an EDCA prioritization mechanism causes a new 
type of coupling between the evolutions of per-class remaining 
idle times and backoff stages, which has been an intricate 
complication lfT9l . another penetration, also formalized by 
mean field argument, has led us to an extended form of an 
ODE model spinning off a generalized FPE as well. 

Lastly, we conjecture that (IMINTb implies the global sta- 
bility of d25l l and d2"6i i as well, as observed in our exhaustive 
simulations. We believe that it is provable with a Lyapunov 
function though the form of which is unknown yet. Although 
theoretical support to this conjecture is not available, we hope 
the discussion can introduce the challenging side of the open 
stability problem. 

References 

[1] M. Benai'm and J.-Y. Le Boudec. A class of mean field interaction 

models for computer and communication systems. Perf. Eval., 65(11- 

12):823-838, Nov. 2008. 
[2] G. Bianchi. Performance analysis of the IEEE 802.11 distributed 

coordination function. IEEE J. Select. Areas Commun., 18(3):535-547, 

Mar. 2000. 

[3] G. Bianchi, I. Tinnirello, and L. Scalia. Understanding 802.1 le 
contention-based prioritization mechanisms and their coexistence with 
legacy 802.11 stations. IEEE Network, 19(4):28-34, July 2005. 

[4] P. Billingsley. Convergence of Probability Measures. Wiley-Interscience, 
2nd ed., 1999. 

[5] C. Bordenave, D. McDonald, and A. Proutiere. Performance of random 
medium access control: An asymptotic approach. In Proc. ACM 
Sigmetrics, June 2008. 

[6] C. Bordenave, D. McDonald, and A. Proutiere. A particle system in 
interaction with a rapidly varying environment: Mean field limits and 
applications. Networks and Heterogeneous Media, 5(l):31-62, Mar. 
2010. 

[7] M. Carvalho and J. Garcia-Luna-Aceves. A scalable model for channel 
access protocols in multihop ad hoc networks. In Proc. ACM MobiCom, 
Sept. 2004. 

[8] K. Duffy. Mean field markov models of wireless local area networks. 

Markov Processes and Related Fields, 16(2):295-328, 2010. 
[9] S. Ethier and T. Kurtz. Characterization and convergence. Wiley, 1986. 



14 



[10] M. Garetto, T. Salonidis, and E. Knightly. Modeling per-flow through- 
put and capturing starvation in CSMA multi-hop wireless networks. 
ACM/IEEE Trans. Networking, 16(4):864-877, Aug. 2008. 

[11] A. Kumar, E. Altaian, D. Miorandi, and M. Goyal. New insights from 
a fixed-point analysis of single cell IEEE 802.11 WLANs. ACM/IEEE 
Trans. Networking, 15(3):588-601, June 2007. 

[12] P. Kumar. An interview with Dr. P. R. Kumar. Science Watch Newsletter, 
available at http:// esi- topics.com/ wireless/interviews/PRKumar. html 
June 2006. 

[13] T. Kurtz. Strong approximation theorems for density dependent Markov 
chains. Stochastic Processes and their Applications, 6:223-240, 1978. 

[14] B.-J. Kwak, N.-O. Song, and L. Miller. Performance analysis of 
exponential backoff. ACM/IEEE Trans. Networking, 13(2):343-355, 
Apr. 2005. 

[15] K. Medepalli and F. Tobagi. Towards performance modeling of IEEE 
802.11 based wireless networks: A unified framework and its applica- 
tions. In Proc. IEEE Infocom, Apr. 2006. 

[16] A. Proutiere. Pushing mean field asymptotics to the limits: Stability 
and performance of random medium access control. In Proc. Stochastic 
Networks Conference, ENS, Paris, June 2008. 

[17] V. Ramaiyan, A. Kumar, and E. Altman. Fixed point analysis of single 
cell IEEE 802. lie WLANs: Uniqueness, multistability. ACM/IEEE 
Trans. Networking, 16(5):1080-1093, Oct. 2008. 

[18] T. Sakurai and H. Vu. MAC access delay of IEEE 802.11 DCF. IEEE 
Trans. Wireless Commun., 6(5):1702-1710, May 2007. 

[19] G. Sharma, A. Ganesh, and P. Key. Performance analysis of contention 
based medium access control protocols. IEEE Trans. Inform. Theory, 
55(4):1665-1681, Apr. 2009. 



Appendix 

A. Alternative Proof of Lemma\l\ 

We show the existence and uniqueness of the equilibrium 
point. Differentiating the right-hand side of (flTT i with respect 
to q, we can see that the following equation determines the 
sign of the derivative. 



k-l 




V K V K >y k +3 
l^k=o l^i=0 7 



Consider a proper subsum, 6 K , which can obtained by replac- 
ing K with k £ {1, • • • , K — 1}. Recall that q > q\ by the 
assumption; then it is easy to see that 8\ < is true. Now 
suppose S K is zero or negative. We show <5 K+1 < if 5 K < 0. 
Rearranging terms of S K+ i, it is not difficult to obtain: 



Of-- — 

' \<ii <?K+1 



where the second term on the right-hand side is zero or 
negative as qi is nonincreasing for i £ {0, • ■ ■ , K}. As S(K) 
is zero or negative, we can conclude that the right-hand side 
of ( fTTT ) is a nonincreasing function which is positive and 

converges to (K + 1)/ EfcLo 9/7 1 at P = 00 • This conclusion 
taken together with the fact that the left-hand side of ( fTTT ) is 
an identical function from [0,oo) to [0,oo) proves that there 
exists a unique equilibrium point q. 

B. Proof of Lemma \3\ 

First, we note from d3TT l that the right-hand sides of (|27| i 
and ( f28l l are nonincreasing in 7 H and 7 C , respectively. The 
proof of this fact is almost identical to that of Lemma Q] 

Assume that there are two solutions (q 11 , q h ) and (q H , q L ) of 
the fixed point equation (T27l)-(|30T> and q H > q H , without loss of 



generality. If we assume that q h > q h , it follows from d30b that 
7 C > 7 ■ Because the right-hand side of ( f28b is nonincreasing 
in 7 C , we must have q L < q h and hence j c > j c . Now we 
have shown by contradiction that g H > q 11 implies q h < q L 
and 7° > 7 C . 

Moreover, we can rewrite < f29b in the following form: 



l-7 R ) A +7 R Ef=o 1 (l 



A-l 



r 



i 



(1-7 R ) £ 



(1-7 R ) A 
7 C 



(36) 



where the second equality can be easily verified. As 7 R = 
1 — is increasing in q H , q n > q 11 implies j c > 7° and 
7 > 7 R - Combining these with the fact that (l36l l is increasing 
in 7 R and 7°, we can establish that g 11 > q H implies 7 11 > 
7 H . On the other hand, since the right-hand side of ( f27T > is 
nonincreasing in j H , the inequality j H > j H must imply g H < 
<i H 

In conclusion, if we assume <7 H > g H , we have q H < q H , 
which implies that q H — q n . Then it automatically follows that 
q L = q L , 7 H = 7 H , and j L — 7 L . 

We have yet to establish the existence of the solution. We 
first note that the left-hand sides of ( f2Tb and (f28T > are identical 
functions of q H and q h , respectively, from [0, 00) to [0,oo). 
Because d30T > is increasing in q h , for each fixed g H , the right- 
hand side of (l28l l is a positive nonincreasing function of q h 
by the proof of Lemma Q] Likewise, as ( f29T > is increasing in 
q H for each fixed q h , the right-hand side of ( |27] i is a positive 
nonincreasing function of <7 H by the proof of Lemma [TJ This 
completes the proof. 

C. Proof of Lemma [4] 

Multiplying both sides of (|27| | and (f28b respectively by (1 — 
7 H ) and (1 — 7°) yields the following equations: 

« H (l-J H ) = * H ^ 0[ l4_ -(1-7 H ), 07) 

Efc=0 „H 



? L (1-7 C ) 



v^^ L t c\ k 
t l E^o (7 C ) 

if L ( 7 c:) fc 
k=0 



(1-7 C )- 



(38) 



The proof is similar to that of Lemma [3] except that: 

(i) We use the fixed point equation d37]l, CH}, (EU and (f30t . 

(ii) We note from Lemma[2]that the right-hand sides of (|37| ) 
and (l38l are decreasing respectively in q 11 and q h , and less than 
or equal to the left-hand sides of ( f37l > and ( f38l > respectively at 
q H = 1 and q L = 1. 

To complete the proof, it is sufficient to show that the left- 
hand sides of d37l i and d38l are increasing respectively in q H 
and q h . It follows from the proof of Lemma|2]that ( f33b implies 
q H < 1 and q h < 1. It is also obvious from the form of 
g L (l - 7 C ) = g L e _<? ""« L that the left-hand side of (O in 
increasing in q h G [0, 1]. 

To sum up again, it is now enough to show that the left- 
hand side of (T37T > is increasing in q H £ [0, 1]. To establish this, 
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we rewrite ( f29T > in a compact form 

V<? H /i(g H ,<f) e~^^ 



7 H = 1 



/l(g H , q L ) + 1 h(q H , q L ) + 1 j 

where h(q H ,q L ) := (e^ A - lj ■ ■ Differentiating 

q H (l — 7 H ) with respect to q 11 yields 



-9 H -9 L U e-« H -e"« H -« L dft 



h + 1 



q 



(h + l) 2 dq H 

where h is a shorthand notation for h(q H ,q h ). The first term 
of the above equation is positive for q H G (0, 1). The sign of 
the second term is determined by which is nonnegative 
because h can be rearranged as 

%- H ,g L ) = ^ C ^.(i-e-^ L ) 
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which is nondecreasing in g H . This completes the proof. 
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